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1.  Introduction 


Developing  sensing  capabilities  for  the  detection  of  ground  targets  concealed  in  a 
forest  environment  necessitates  a  fundamental  understanding  of  the  scattering 
properties  of  not  only  the  targets  and  the  trees  but  also  their  interactions  over  the 
observation  time  interval.  For  a  given  target  in  motion  under  the  tree  canopy,  the 
time-frequency  behaviors  of  these  scattering  interactions  are  very  much  dependent 
upon  the  tree  arrangement  and  the  precise  structure  of  the  canopy.  In  this  respect, 
the  spatial  distribution  of  the  various  tree  constituents  (trunks,  branches,  and 
foliage)  potentially  becomes  an  important  determinant  of  the  detection  performance 
of  the  radar  sensor  and  processing  algorithm,  and  therefore,  needs  to  be  preserved 
in  a  wave  scattering  solver.  For  example,  in  the  application  of  a  clutter-limited 
technique  such  as  along-track  interferometry  for  moving  target  indication,  the  trees 
surrounding  the  target  can  introduce  deleterious  effects  into  the  interferometric 
phase,  leading  to  decorrelation  and  a  nonlinear  target  phase-velocity  relationship, 
as  a  result  of  tree  occlusion,  attenuation,  and  phase  distortion.  From  the  radar 
simulation  perspective,  such  effects  can  only  be  studied  with  a  physics-based  solver 
that  faithfully  captures  the  actual  positions  of  the  target  and  the  trees  and  their 
scattering  phase  centers  over  the  entire  coherent  processing  interval. 

Existing  foliage -penetration  radar  performance  analysis  tools  are  rather  simplistic 
in  that  they  do  not  take  into  account  the  localized  interactions  of  the  target  with  its 
surrounding,  nor  do  they  fully  consider  the  polarimetric  scattering  responses  of  the 
scene  (i.e.,  the  target  and  clutter  responses  are  often  rudimentarily  incorporated  into 
the  signal  model  just  by  using  radar  cross  section  values  or  coefficients).  Higher 
fidelity  scattering  analysis  in  the  form  of  the  discrete  scatterer  technique  (DST)  has 
been  undertaken  in  the  microwave  remote  sensing  community  within  the  context 
of  developing  algorithms  for  retrieving  tree/forest  biophysical  parameters  from  the 
scattering  returns.1,2  In  Lin  and  Sarabandi,1  wave  theory  is  applied  in  constructing 
a  3-D  coherent  scattering  model  for  forest  stands  consisting  of  fractal-generated 
trees;  unfortunately,  in  this  model,  although  individual  trees  have  deterministic 
positions,  the  propagation  or  transmissivity  effects  of  the  trees  are  considered  only 
in  a  statistical  manner  within  a  layered-canopy  representation.  A  similar  approach 
to  that  of  Lin  and  Sarabandi1  is  also  employed  by  Thirion  et  al.,2  wherein  the  main 
physical  features  of  the  canopy  such  as  the  trunks  and  crown  shapes  are  captured, 
but  the  branches  and  leaves  are  specified  only  indirectly  through  probability 
distributions. 

Moreover,  notwithstanding  the  usefulness  of  the  analytical  solvers  developed  in 
these  aforementioned  studies,  and  others  not  discussed  here,3-6  the  application 
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focus  in  existing  works  has  not  been  on  the  detection  of  moving  targets  under  the 
canopy  but  rather  only  on  the  responses  of  the  trees  themselves.  It  should  be  noted 
that  a  variety  of  radiative  transfer-based  solvers  also  have  been  put  forth  for  tree 
response  characterization.  While  traditional  radiative  transfer  theory  has  been 
applied  successfully  to  a  number  of  remote  sensing  scenarios — some  involving 
complex  forest  canopy  structures7 — this  method  does  not  include  the  coherent 
scattering  effects  occurring  within  the  tree  structure  and  therefore  cannot  be  used 
to  compute  the  phase  response  of  a  scene. 

In  order  to  account  for  the  spatial  distribution  of  the  trees  coherently  and 
systematically,  this  work  follows  the  framework  proposed  by  Liu  et  al.,8  in  which 
a  multi-ray  approach  is  exploited  to  calculate  the  aggregate  scattered  fields  from  a 
collection  of  trees,  with  the  transmissivity  of  each  ray  determined  from  a  cell-based 
representation  of  the  scene.  A  cell-by-cell,  interaction-tracking  procedure  is  also 
proposed  here  for  speeding  up  the  calculations  of  the  responses  from  small 
components  such  stems  and  leaves,  which  are  expected  to  far  outnumber  those 
scatterers  composing  the  main  trunks.  Furthermore,  tree  and  scene  generation  for 
the  analytical  solver  is  carried  out  in  conjunction  with  that  for  a  large-scale,  full- 
wave  solver,  allowing  for  an  assessment  of  the  accuracy  of  the  ray-based  solution 
for  any  scene.  Note  that  the  details  of  the  full- wave  solver  itself  have  been  published 
elsewhere,9  which  demonstrated  the  practicality  of  finite-difference  time-domain 
(FDTD)-based  scattering  simulations  of  a  realistic  forest  scene  for  foliage- 
penetration  radar  applications. 

In  brief,  this  work  is  organized  as  follows.  In  Section  2,  an  overview  of  the 
procedure  for  generating  3-D  realistic  trees  for  both  the  analytical  and  numerical 
solvers  is  presented.  In  Section  3,  the  formulations  for  the  DST  are  described  in 
detail.  In  Section  4,  computational  examples  are  given:  the  DST  solutions  are 
compared  with  full-wave  simulation  results;  the  accuracy  of  the  cell-by-cell 
approach  is  evaluated;  and  a  general  discussion  on  the  scattering  effects  of  a  forest 
stand  is  put  forth.  Finally,  Section  5  summarizes  the  study. 

2.  Tree  Generation 


Realistic  3-D  tree  structures  are  created  using  the  open-source  tree  generation 
engine  Arbaro,10  which  implements  the  recursive,  rule-based,  growth  algorithm 
developed  by  Weber  and  Penn.1 1  Given  a  random  seed  and  an  input  parameter  file 
defining  the  physical  characteristics  of  a  tree  (e.g.,  tree  and  trunk  shape,  number  of 
branch  levels,  branch  sizes  and  angles,  leaf  shape  and  density),  the  algorithm 
constructs  the  hierarchical  structure  for  the  tree  by  following  a  series  of  recursive 
rules  derived  from  geometrical  observables.  For  the  DST  solver  of  this  work,  each 
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tree  is  exported  from  Arbaro  as  a  multilevel  Persistence  of  Vision  Raytracer 
(POV-ray)  file  in  which  the  trunk/branches  and  leaves  are  represented  by  tube  and 
disk  primitives,  respectively.  Additional  geometric  transformation  or  processing  of 
the  exported  tree  object  (i.e.,  decomposition  of  the  tube  primitives  into  cylinders  of 
desired  sizes,  global  translation,  and  scaling,  etc.)  is  carried  out  before  the  entire 
model  is  passed  on  to  the  scattering  solver.  For  the  FDTD-based  simulations,  the 
mesh  form  of  the  tree  structure  is  first  exported  from  Arbaro  as  a  Wavefront  OBJ 
file;  subsequently,  voxelization  of  the  mesh  is  performed  using  the  parity  count  and 
ray  stabbing  methods. 9,12-13 


3.  Cell-Based  Discrete  Scatterer  Technique 


Once  the  trees  have  been  parsed  into  canonical  components  or  scatterers  (i.e., 
dielectric  cylinders  and  disks),  each  with  deterministic  locations  and  well-defined 
physical  and  electrical  properties,  the  scattering  solver  is  derived  by  applying  the 
3-step  procedure  as  follows.  First,  the  tree  scene — consisting  Nt  trees — is 
discretized  into  Nt  cubic  cells  (Fig.  1),  with  each  scatterer  assigned  to  a  cell. 
Second,  the  propagation  paths  for  each  scatterer  are  identified  and  tabulated:  this 
necessitates  deducing  all  the  cells  intercepted  by  the  4  rays  shown  in  Fig.  2.  Finally, 
the  ray-based  scattering  response  for  each  scatterer  is  calculated  analytically,  while 
taking  into  consideration  the  effects  of  attenuation  and  phase  change  encountered 
along  each  ray  path,  and  then  the  total  scattering  response  of  the  entire  scene  is 
obtained  by  coherently  summing  the  returns  from  all  the  scatterers.  Succinctly,  in 
the  multi-ray  approach,  for  incident  direction  ki  and  scattering  direction  ks ,  the 

scattering  matrix  for  the  nth  («  G  {l,2 ,...,NC})  scatterer  can  be  written  as 


where 


(*,);  (2) 


leg 


(3) 

(4) 
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9  gcg 


■Te 


(5) 

in  which  p"  is  the  location  of  the  scatterer  center;  r"  -  r"  -  2 (ro"  ■  zjz ; 

kgi  =  ki  -  2\ki  •  z  jz ;  kgs  =  ks  -  2 \ks  -zjz;  S  is  the  scatterer  free-space  response; 
— « 

Tu  (•)  (n  =  c,  g,  eg,  gc)  is  the  transmissivity  matrix  of  the  propagation  path  as 
determined  by  the  cells  intercepted;  and  r(-)  is  the  ground  reflection  coefficient 


matrix.  Essentially,  in  the  presence  of  the  ground,  the  multi-ray  technique  captures 
the  first  4  orders  of  scattering  effects:  Eq.  2  is  the  direct  scattering  response  from 
the  scatterer;  Eqs.  3  and  4  correspond  to  the  single  ground-bounce  returns;  and 
Eq.  5  is  the  double  ground-bounce  contribution. 


Fig.  1  Cellularized  DST  computational  domain 
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(a) 


(b) 


k, 


(c) 


(d) 


Fig.  2  Ray  paths  for  a  scatterer:  a)  direct  response,  b-c)  single  ground-bounce  returns,  and 
d)  double  ground-bounce  return 

=  fs,n 

In  Eqs.  2-5,  the  free-space  response  of  the  scatterer,  S  (•,•),  is  derived  with  the 
application  of  2  sets  of  approximate  fonnulations:  for  a  finite-length,  circular 
cylinder  scatterer,  the  volumetric  current  integration  method  is  employed  wherein 
the  internal  cylinder  fields  are  assumed  to  be  the  same  as  those  of  the  infinite-length 
case;14,15  and  for  a  thin-disk  scatterer,  the  generalized  Rayleigh-Gans  solution  is 
applied  wherein  the  internal  disk  fields  are  related  to  the  incident  fields  through  the 
polarizability  tensor  of  the  scatterer.1’’  In  the  current  work,  the  leaves  are  simply 
modeled  as  circular  disks.  The  analytical  formulations  for  these  2  sets  of  solutions 
are  included  in  the  Appendix. 

The  aforementioned  transmissivity  matrix 


(6) 


is  calculated  with  Foldy’s  approximation,16  which  derives  the  distortion  effects  of 
a  random  medium  on  a  mean  field.  It  can  be  shown  that 


and,  assuming  a  medium  with  azimuthal  symmetry,  T"hv(-)=  7^v/i(-)  =  0.  In  Eq.  7, 
'b  is  the  propagation  path  length  vector  (i.e.,  the  bth  element  of  !"(•)  is 
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the  path  length  of  the  ray  in  cell  b,  b  £  {l,2,...,Nb});  and  Kpp (•)  e  C  v'’  is  the 

aggregate  forward  scattering  vector  for  all  the  cells,  with  its  hth  element  computed 
as 


40 =j~  Z  4’ (•■•).  (8) 

where  v  is  the  cell  volume.  Together,  Eqs.  1-8  give  the  most  general  framework  for 
analyzing  the  scattering  from  a  tree  scene,  taking  into  account  propagation 
attenuation  and  phase  change  effects  experienced  by  the  incident  and  scattered 
waves  for  each  discrete  scatterer  in  the  random  medium.  In  view  of  the  fact  that  the 
scatterers  that  compose  the  tree  trunks  are  the  largest  scatterers  in  the  scene  and 
therefore  would  provide  the  dominant  responses,  extra  consideration  must  be  taken 
in  evaluating  Eq.  6  for  the  trunk  elements.  In  most  realistic  scenarios,  the  individual 
scatterers  of  a  trunk  are  not  expected  to  be  shadowed  by  other  scatterers  belonging 
to  the  same  trunk — in  other  words,  a  trunk  cannot  shadow  itself.  As  such,  in 
calculating  the  transmissivity  matrices  for  the  scatterers  of  a  particular  tree  trunk  t 
(«etmnkt,fe{l,2,...,A,}),  the  effects  of  self-shadowing  should  be  excluded  by 

modifying  Eq.  8  as 

40=  Z  C (•■•)•  (?) 

mecell  b 

m  g  trunk  t 

To  increase  the  computational  efficiency  of  the  solver,  instead  of  a  scatterer-by- 
scatterer  brute  force  approach,  a  cell-by-cell  scattering  and  propagation  strategy  can 
be  applied  for  the  treatment  of  the  smallest  branch  elements  and  the  leaves,  which 
are  expected  to  far  outnumber  those  scatterers  composing  the  trunks  and  larger 
branches.  In  this  approach,  all  the  constituents  of  a  cell  are  assumed  to  have  the 
same  set  of  transmissivity  matrices  (for  each  ray  path)  and  the  responses  are 
propagated  together  as  a  group  rather  than  one  by  one.  To  this  end,  the 
approximation  is  made  such  that  the  ray  paths  of  all  the  scatterers  belonging  to  the 
same  cell  are  identical.  Here,  the  center  of  the  cell  is  used  as  the  common  reference 
point  for  identifying  all  those  paths  and  evaluating  the  transmissivity  matrices  using 
Eqs.  6-8.  The  accuracy  of  the  previous  approximation  is  investigated  in  the  next 
section. 

Essentially,  in  this  cell-by-cell  strategy,  computational  speedup  is  provided  by  not 
having  to  calculate  the  transmissivity  matrices  on  a  scatterer-by-scatterer  basis,  as 
well  as  by  not  having  to  determine  and  tabulate  the  ray  paths  for  every  scatterer  in 
a  cell.  Note  that  the  same  cell-based  strategy  can  potentially  be  applied  to  the 
scatterers  for  the  trunks  and  larger  branches,  as  well.  However,  in  general,  since 
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those  scatterers  are  expected  to  provide  the  more  dominant  responses — as 
compared  to  the  smallest  branches  and  the  leaves,  the  associated  errors  induced  by 
the  underlying  assumptions  of  the  technique  are  also  expected  to  be  more 
significant. 

4.  Numerical  Results 


The  scenes  of  interest  are  displayed  in  Fig.  3,  which  shows  a  forest  stand  composed 
of  Nt  =  72  randomly  generated  quaking  aspen  ( Populus  tremuloides )  trees  modeled 
with  different  levels  of  structural  fidelity.  The  trees  are  constructed  assuming  a 
4-level  hierarchy — as  such,  there  are  4  cases  of  consideration:  the  trunks  only,  the 
trunks  with  primary  branches,  the  trunks  with  primary  and  secondary  branches,  and 
the  trees  with  all  their  branches  and  leaves.  As  the  starting  point,  the  backscattering 
returns  of  the  forest  stand  are  first  examined  by  comparing  the  DST  results  with 
those  from  the  full-wave  solver.  Note  that  since  the  FDTD  method  as  employed 
here  cannot  efficiently  model  extremely  thin  structures  such  as  leaves,  validation 
of  the  DST  results  is  only  carried  out  for  scene  configurations  without  the  foliage. 
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Fig.  3  Forest  stand  with  72  randomly  generated  trees:  a)  tree  trunks  only,  Nc  =  2017,  b) 
tree  trunks  +  primary  branches,  Nc  =  25397,  c)  tree  trunks  +  primary  and  secondary  branches, 
Nc  =  292877,  and  d)  tree  trunks  +  primary  and  secondary  branches  +  leaves,  Nc  =  1878880 


The  FDTD-  and  DST-calculated  backscattering  radar  cross  sections  o- 
(p,q  =  v,h)  and  phase  differences  (d  =  c,x)  of  the  forest  stand  are  illustrated  in 

Figs.  4-6,  as  a  function  of  the  elevation  incidence  angle  9i  and  structural  fidelity. 
Formally,  the  plotted  co-  and  cross-polarized  responses  are  defined  as  follows: 


°M=4;r 


pq 


£  =  arg  (sji, 


,wLJhh  /  ’ 


and  &  =  arg (SVVS, 


,vvLJ hv  /  ? 


where 


is  the 


/^/-polarized  scattering  matrix  element.  Note  that  in  the  backscattering  direction, 
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Shv  -  ~Svh .  The  results  at  each  elevation  incidence  angle  are  obtained  by  averaging 

the  responses  over  36  evenly  spaced  azimuth  incidence  angles  spanning 
(f)L  =  [0°,  360°]  and  over  301  frequencies  spanning  /  =  [200  MHz,  500  MHz] 
(P-band).  The  ground  and  the  tree  trunk/branches  are  assumed  to  be  electrically 
homogeneous,  with  relative  dielectric  constant  and  conductivity  ( sr ,  a d)  of 
(5.45,  20  mS/m)  and  (13.9,  39  mS/m),  respectively.  The  FDTD  computational 
domain  has  dimensions  30  m  x  30  m  x  9.2  m  and  is  discretized  into  2  billion  voxels; 
parallelized  simulations  are  performed  after  partitioning  each  scene  into  256 
subdomains.9  On  the  other  hand,  the  DST  scene  has  dimensions  30  m  x  30  m  x 
10  m  and  is  divided  into  Nb  =  9,000  cubic  cells  (i.e.,  each  cell  has  a  side  length  of 
1  m). 


Fig.  4  Comparison  of  DST  solution  with  FDTD  solution  for  scene  containing  trunks  only 
(Fig.  3a):  a)  radar  cross  sections  and  b)  phase  differences 
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Fig.  5  Comparison  of  DST  solution  with  FDTD  solution  for  scene  containing  trunks  + 
primary  branches  (Fig.  3b):  a)  radar  cross  sections  and  b)  phase  differences 
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Fig.  6  Comparison  of  DST  solution  with  FDTD  solution  for  scene  containing  trunks  + 
primary  and  secondary  branches  (Fig.  3c):  a)  radar  cross  sections  and  b)  phase  differences 

As  evident  from  Figs.  4-6,  in  general,  except  for  the  cross-polarized  response  of 
the  trunks-only  scene  (Fig.  4),  the  DST  results  are  in  good  agreement  with  those 
from  the  full-wave  method — in  both  amplitude  and  phase.  For  the  trunks-only 
scene,  the  observed  discrepancy  between  the  analytical  and  numerical  cross- 
polarized  returns  is  in  part  due  to  the  fact  that  the  FDTD  solver  as  implemented 
here  cannot  capture  low  cross-polarization  signal  levels.  Once  the  branches  are 
introduced  into  the  scene,  the  cross-polarized  response  increases;  therefore,  those 
cases  with  the  branches  (Figs.  5  and  6)  do  not  pose  a  problem  for  the  FDTD  solver. 
Note  that  at  the  cost  of  increased  computational  time  and  memory,  the  fidelity  of 
the  numerical  cross-polarized  solution  can  be  improved  by  reducing  the  FDTD 
voxel  size.  Although  the  limitation  of  the  full-wave  analysis  here  precludes  a 
complete  validation  of  the  DST  results  when  the  cross-polarized  response  is  low, 
the  DST  co-polarized  solutions  are  seen  to  match  those  from  FDTD  reasonably  well 
in  all  cases. 
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The  analytical  solutions  presented  thus  far  are  obtained  by  calculating  the  responses 
from  the  scene  on  a  scatterer-by-scatterer  basis.  To  speed  up  the  computations,  the 
cell-by-cell  interaction  algorithm  discussed  in  Section  3  is  applied  for  the  treatment 
of  the  secondary  branches  in  the  scene  in  Fig.  3c.  A  comparison  of  the  results  from 
this  new  approach  with  those  from  the  original  scatterer-by-scatterer  technique  is 
shown  in  Fig.  7.  It  is  observed  that  the  more  computationally  efficient  approach 
retains  the  accuracy  of  the  brute  force  technique.  Note  that  the  memory  usage  and 
runtime  for  the  scatterer-by-scatterer  simulation  are  369  MB  and  33  min/frequency, 
respectively,  on  a  Dell  Precision  T7500  workstation  with  Intel  Xeon  CPU  of 
2.67  GHz,  while  those  for  the  cell-by-cell  simulation  are  200  MB  and 
40  s/frequency,  respectively. 
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-30 

'350  15  30  45  60  75  90 
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Fig.  7  Comparison  of  scatterer-by-scatterer  vs.  cell-by-cell  DST  solutions  for  scene 
containing  trunks  +  primary  and  secondary  branches  (Fig.  3c):  a)  radar  cross  sections  and 
b)  phase  differences 
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The  DST-simulated  backscattering  responses  of  the  forest  stand  for  the  4  levels  of 
scene  structural  fidelity  are  plotted  together  in  Figs.  8  and  9.  The  foliage  is  modeled 
by  considering  the  leaves  as  circular  dielectric  disks  with  a  thickness  of  0.2  mm,  a 
mean  radius  of  2.6  cm,  and  (sr,  ad)  =  (28.0,  0.37  S/m).  The  cell-by-cell  strategy  is 
implemented  for  the  secondary  branches  and  the  leaves  when  the  scene  contains 
those  components.  (The  memory  usage  and  runtime  for  the  case  with  all  the 
branches  and  leaves  are  1.0  GB  and  56  s/frequency,  respectively.)  For  the 
co-polarized  responses,  it  is  noted  that  as  the  elevation  incidence  angle  increases, 
the  dominant  contribution  tends  to  be  provided  by  the  tree  trunks — though  this  is 
more  apparent  for  hh  than  vv.  In  general,  the  inclusion  of  the  primary  branches  is 
seen  to  lead  to  an  enhanced  backscattering  intensity,  compared  to  the  response  of 
the  scene  with  only  the  trunks.  However,  further  addition  of  the  secondary  branches 
and  the  leaves  to  the  scene  can  weaken  the  return — especially  at  the  higher 
incidence  angles  where  the  propagation  loss  of  the  canopy  plays  a  more  important 
role.  Unlike  the  co-polarized  returns,  the  cross-polarized  response  of  the  forest 
stand  is  mainly  determined  by  the  primary  branches  and  not  by  the  trunks.  This  is 
as  expected,  since  the  trunks  here  resemble  slightly  asymmetric  cones  that  have 
only  a  small  cross-polarized  backscattering  response.  As  in  the  co-polarized  cases, 
the  inclusion  of  the  secondary  branches  and  the  leaves  tends  to  decrease  the  cross- 
polarized  return,  relative  to  that  for  the  scene  with  only  the  trunks  and  the  primary 
branches. 
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Fig.  8  Radar  cross  sections  of  the  scenes  in  Fig.  3:  a)  ou,  b)  a i,i„  and  c)  ai,v 
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Fig.  9  Phase  differences  of  the  scenes  in  Fig.  3:  a)  £  and  b)  £v. 

5.  Conclusions 


An  analytical  solver  is  developed  in  this  study  for  characterizing  the  far-field 
electromagnetic  scattering  responses  of  trees.  The  solution  is  based  on  a  multi-ray 
discrete  scatterer  approach,  applicable  to  a  cellularized  scene  comprising  realistic 
3-D  tree  structures  of  which  the  constituents  are  represented  by  cylinders  and  disks. 
The  accuracy  of  the  solver  is  determined  by  comparing  the  backscattering  results 
for  a  forest  stand  at  P-band  with  those  from  large-scale,  full-wave  simulations.  A 
cell-by-cell  strategy  is  also  proposed  and  implemented  to  increase  the 
computational  efficiency  of  the  solver  for  evaluating  the  scattering  contributions 
from  the  smallest  branches  and  the  leaves.  Overall,  the  analysis  and  computational 
framework  put  forth  is  intended  as  a  tool  for  investigating  scattering  effects  relevant 
to  airborne  foliage-penetration  radar  applications. 
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Appendix.  Scattering  Solution  for  Branch  and  Leaf  Scatterers 


Approved  for  public  release;  distribution  is  unlimited. 


19 


Analytical  fonnulations  for  calculating  the  electromagnetic  scattering  responses  of 
tree  branches  and  leaves  are  given.  Note  that  the  following  treatment  assumes  the 
incident  and  scattered  wave  directions  ( k: ,  ks )  and  their  associated  polarization 
vectors  (  h  ,  v)  are  defined  by 

ki  =  sin  0!  cos  (jix  +  sin  0j  sin  fay  -  cos  0  z\  (A- 1 ) 

ks  =  sin  0s  cos  (j)sx  +  sin  9S  sin  <f>sy  +  cos  0z\  (A-2) 


i  _ 

Us  |  zxk. 
v.  =  h  x  k  ; 

l,S  l.S  l,S> 


(A-3) 

(A-4) 


in  which  the  subscripts  i  and  s  identify  quantities  related  to  the  incident  and 
scattered  waves,  respectively;  6i  is  the  incidence  angle  in  elevation  (measured  from 
-  z  );  Os  is  the  scattering  angle  in  elevation  (measured  from  +  z  );  and  <fc,s  is  the 
incidence/scattering  angle  in  azimuth  (measured  from  +  x  ). 


A-l  Scattering  from  Branches 


In  this  study,  the  tree  branches  are  approximated  with  circular  cylinders.  The  free- 
space  scattering  matrix  elements  of  a  finite-length,  dielectric  circular  cylinder  with 
length  L,  radius  a,  and  complex  relative  dielectric  constant  er,t  are  derived  as1,2 


,((,(■)  = 


Ps  ■- 


kl  (er  l  —  l)  L  (  kaL  ( cos  Ot  +  cos  9S ) 


+o°  _ 


K , 


jn{fs~fi) 


■>  p?q  =  v,h 

(A- 5) 


where 


k,,„  =  /kJ2f' '  (A,  -jccn){x-jy)-  k2,„=-  2k 


(A-6) 


For  q  =  v, 


iM  sin  0  cos 0  n  nW sin  0.  cos 0  jM  sin 0 

jy  —  J  n _ i _ l_  .  n  —  _ i _ l_  .  y  _  J  n  l 

7  x  ’  rn  T-  (  i  \  ^  9  in 


J„(kp,a)R„ 


Jn(kp,ra)Rn  ’  "  Jn(kp^)Rn  ’ 


(A-7) 
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for  q  =  h, 


n  W  sin  0,  cos2  0,  a  ~.iN„  sin  0  nW sin  <9.  cos  <9. 

= - 7 - 7 - ;  Pn  =  - r— ;  Yn  = - ' 


Jn  (kp,ra)R„ 


Jn(kP,ra)Rn  ’  "  Jn(kpra)Rn 


(A-8) 


and 
/  =- 


(v°r  ~{kp,sa) 


\kp,raJ„  [kPpa)J„+l  {kp,ra)-kp,SaJ„  (kp,ra)Jn+ 1  (**,«)); 


(A-9) 


7r(knia)  Hi2)  (k  a)  r 
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M  = 


N  = 


gjAM  -t. ’(*,„*) 


kn  aH{2)  [k„  a)  A  a/  (A„  a) 

n  y  p,i  J  p,r  n  \  p,r  ) 

n? J.  iK.fi) 


kf.fiH',:'{kp.fi)  "  kP.,aJ,  (kA 

1  1 


W  = 


(kp,ia)2  {kp,raf 


(A-l  1) 

(A- 12) 
(A- 13) 


The  previous  equations  assume  kp  .  =  k0  sin  o, ,  kps=ko  sin  0s ,  and 

kp,r=ko^£r,,~COs20,  ■ 

A-2  Scattering  from  Leaves 

The  scattering  response  of  a  circular  disk-shaped  leaf  with  thickness  r,  radius  a,  and 
complex  relative  dielectric  constant  e,-.t  is  given  by 

spq  [k  A)  =  —  J^S,a\Sr ,  -  l)  [aT  {PsAi)  +  iaN  -  aT  ){Ps  ■ z)  (*  •  Vi  )]  »  PA  =  A  h 

U  Hsia 


(A- 14) 

where 

am 

Vr  =  ; 

c  2 

(A- 15) 

1 

“T  ~  (su-l)gI+l’ 

(A- 16) 
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(A- 17) 


(A- 18) 

(A- 19) 
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1 


2  a 

m  =  — , 
T 


Qsi  =  ko  (sin2  0S  +  sin2  6>  -  2 sin  0s  sin  6>.  cos(A  -  r/t  ))-  . 
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